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(54) Method of performing geometric measurements on digital radiological images using 
graphical templates 

(57) Method of performing geometric measure- 
ments on digital radiological images, using graphical 
templates. . . 

A digital image representation of a radiological im- 
age is acquired and displayed on a display device. A 



it template, comprising the set of measure- 
ment points or areas to be measured, is retrieved and 
displayed adjacent to the digital image. The. measure- 
ment points or objects are mapped on to the displayed 
image and the resultant measurements stored; 
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Description 

FIELD OF THE INVENTION _ 

[0001] The present invention relates to computer-assisted radiological measurements on radiographic images. 

BACKGROUND OF THE INVENTION 

[0002] In musculoskeletal radiology, it is current practice to use a wide range of. imaging modalities to determine 
skeletal disorders and abnormalities. 

[0003] In this discipline diagnosis is often based on quantified radiological findings of geometrical quantities. 

At least 50% of all radiological examinations today-are conventional exposures of thorax and the skeleton. 80% of 

skeletal exposures leads to the correct diagnosis based on the radiographs. 

A plurality of textbooks in the X-ray diagnostics from world-renowned radiological and orthopaedic experts make a 
contribution in this area. 

[0004] However, to refine the diagnosis, to conduct better differential diagnosis, to assess severity of change, to plan 
and control therapy, to conduct treatment follow-up, to establish physical ability in sports medicine, labour medicine 
and military medicine, verbal descriptions based on skeletal radiographs are insufficient in many respects. 
Better diagnosis can be achieved by quantifying radiological findings. 

[0005] Geometrical quantities measured on radiological images must be checked against normal values. These 
normal values have been collected from measurements of a representative sample of the normal healthy population 
and are tabulated in the above-mentioned text books. 

[0006} . Geometrical measurements in digital images comprise linear and angular measurements. Linear measure- 
ments in 2 dimensions and 3 dimensions may be supplemented with distance along a curvilinear path. Angular meas- 
urements are considered in the plane of the image, in a world plane or in 3D space. Geometrical areas are considered 
in the image plane t or more generally of surface patches in 3D images. Volumes are computed in 3D images but may 
be based on planar measurements. Geometrical indices are clinical quantities basedon image measurements. In any 
of these categories, the measurand Is defined as the physical parameter being quantified by measurement. 
[0007] Today, radiological measurements on X-ray images are either made on film using conventional measuring 
devices (such as a ruler, a caliper or a rubber band to measure lengths, and a square orgoniometerto measure angles) 
or in a digital image displayed on screen using cursor controlled points (such as a pair of points to measure Euclidean 
distance between). 

[0008] The current measurement procedure thus involves 4 distinct media: 

1 . An X-ray film comprising the anatomical sites to be measured, displayed on a light box. With the emergence of 
digital radiography modalities (film digitisation, computed radiography/digital radiography sensors) , the digital im- 
age may be displayed on a computer display. However, such electronic medium still is physically different from 
the other components described hereafter. 

2. A measurement atlas comprising the measurement scheme: imaging technique, graphical template and de- 
scription of the measurements covered by the scheme (nomenclature, clinical significance, and normative tables, 

. sometimes interchangeably represented by curves) 

3. An analogue measurement device (ruler, square...) to perform geometrical measurements, 

4. Pencil/Paper to note the measurement quantity according to the appropriate medical nomenclature and.the 
measurement value. 

[0009] A calculator device may be needed to compute indices from a collection of measurements, or to convert 
measured values to true quantities using calibration measures. Alternatively, electronic spreadsheets may be used in 
conjunction with a database to store the "measurements and indices. 

. [001 0] .The use of different media asks for repeated focusing of attention between the atlas and the radiological image. 
Moreover, in the absence of an atlas scheme the position of the measurement objects is riot defined and hence different 
users may locate a given anatomical landmark differently. 

[0011] Because there is no link between the measurement template in the atlas and its associated measurement 
entities, there fe no systematically imposed consistency of the naming of measuredquantities, and therefore exchange 
or collection of measurement values of different clinicians (e.g. for the purpose of cross*refereeing) is fundamentally 
hampered. ' 

[0012] Another major drawback of the prior art method to perform geometricat measurements-is increased meas-" 
urement error or measurement uncertainty. 

The error of measurement is the result of a measurement value minus the (true) value of the measurand. Measurement 
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error is due to different sources, basically falling into one of two classes: systematic and random" ewers.-.- 
[0013] Systematic or bias errors arise from consistent and repeatable sources of error (like an offset in calibration). 
Systematic errors can be studied through inter-comparisons, calibrations, and error propagation from estimated sys- 
tematic uncertainties in the sensors used. Systematic error is defined as the mean that would result from an infinite 
number of measurements of the same measurand carried out under repeatability conditions minus the (true) value of 
the measurand. This source of error can be reduced by better equipment and by calibration. 

[0014] Random errors also referred to as statistical errors, arise from random fluctuations in the measurements. In 
particular, digitisation noise (e.g. geometric digitisation: finite pixel size; intensity digitisation: quantisation of grey levels) 
and the errors introduced by counting finite number of events (e.g. X-ray photon count) are examples of random errors 
in the context of digital X-ray images. Random error is defined as the result of a measurement minus the measurement 
that would result from an infinite number of measurements of the same measurand carried out under repeatability 
conditions. Particularly this source of error is prevailing in the prior art of performing measurements on X-ray images: 
Inter-observer and intra-observer variance on measurement values contribute to this source of error, and has its origin 
in several forms of ambiguity in defining the measurand. ' 
Lack of unambiguous definition of the measurand with respect to the imaged patient anatomy and tack of knowledge 
of the geometrical pose of the patient with respect to source and detector are the main source of random error. 
[0015] Repeatability and reproducibility of a measurement require that the random errors involved in the measure- 
ment procedure are low. Although random errors are reduced when a measurement is repeated many times and the 
results averaged together, this can rarely be achieved in clinical praclice. 

[001 6] It Is an object of the present invention to provide a user-friendly radiological measurement method that over- 
' comes the drawbacks of the prior art. 

SUMMARY OF THE INVENTION 

[0017J The above-mentioned objects are achieved by amethod as set out in claim 1. 

(0018] The method of the present invention is described with regard to geometrical measurements performed or* 
radiological images of humans. It will be clear that the invention is not limited to human beings and can also be-appiied 
in other fields, for example in veterinary applications. 

[0019] In the context of the present application the term 'activation' refers to loading a measurement scheme from 
memory and measurements to be performed according to the loaded scheme. 

[0020J A measurement scheme or measurement template is a pattern of measurements to be performed. The meas- 
urements to be performed are grouped in the form of a measurement procedure wherein the sequence, th&: inter- 
dependence and method of measurements are defined. Such a measurement scheme can be noted and stored in a 
computer in standard notation XML (Extensible Mark-up Language). 

[0021 J In general the measurement scheme comprises a graphical part (also called graphical model) and an internal 
part (also called interna) model). 

[0022] The graphical part represents the geometric relation between measurement entities (objects and operators) 
arid the anatomy in the type of image on which the measurements are to be performed. Measurement objects are e. 
g. points, lines, circles etc. The measurement objects are defined relative to the anatomy. The intended position with 
respect to the anatomy is thereby unambiguously denoted. 

• The measurement objects carrbe labelled e.g. with tfie appropriate medical nomenclature.Otfier kinds of naming are 
possible, e.g. for naming intermediate objects needed in the course of the measurement procedure. : \ 
Measurement-entities also comprise distances, angles which are the result of the measurement procedure and. which 
^^^e^^jriso-n^nedaccorclins to medjcaf,csnyentionv. .. '-'A 

• [0023}- The graphical part may also comprise measurement mark up such as a graphical: indication of distance or 
angle. . : . : ...:... 

It may also comprise textual annotations of the names of objects and entities. . : 

[0024} In the graphical part, a region of interest box around a set of measurement points may be superimposed, to 
indicate the area to which geometric operators (such as zooming) or intensity operators (such as enhancement or 
landmark extraction operators) can be applied. 

To ensure correct matchwith the image, e.g. with the. .imaged patient body part, the correct exposure parameters can : 
•be-specifled In the template, quantifying the imaging geometry. Angle of incidence and source-patient-detector disv < 
tances are the most important geometric parameters, which much be observed at exposure time fora specific template 
to be applicable. Specific data structures and methods of each measurement entity and object class are devoted to 
control their graphical content and behaviour. 

[0025} The measurement scheme may also comprise normative values associated with the mcasuremenientities. ' 
These normative values are.e.g. organised as a normative value table and are used for comparison with measured. :.: 
entities. Ataomiaf values may be signalled. 
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[0026] The internal part represents the functional dependencies between measurement entitteF.-Panctienal depend- 
encies comprise type of measurements, measurement methods, location of a measurement entity in a measurement 
dependency graph, order in which measurements are to be performed etc. _ 
[0027] The internal model can be implemented as an object-oriented model of a coherent set of measurements to 
s be performed on a digital (medical) image. 

It can for example be represented by a measurement dependency graph. The internal model decomposes eachmeas- 
urement of a physical quantity into its constituent measurement objects. It further specifies methods to map the objects 
geometrically onto the medical image. 

In informatics terms, the model is composed of objects belonging to an associated class, consisting of a collection of 
'0 data structures and methods, said methods operating on the corresponding data structures. 

The flow of execution of the measurements of the measurement scheme is imposed by a measurement dependency 
graph, in which the nodes correspond to measurement of geometric objects, and directed arcs define the relationship 
between the nodes. ( 
[0028] As will be explained further on, the internal part may activate the nodes in several ways. 

15 

• sequential point-only operation willfirst address admeasurement points prior to evaluating the measurement object 
and operator nodes; 

• the sequential measurement object operation will let fire each node immediately when the results of all children 
are available, representing the pure data-flow case; 

20 . a node representing a set of points linked to an active contour model will invoke an automatic segmentation com- 
putation in order to compute the location of the image mapped points prior to firing the nodes representing geo- 
metrical objects based on said mapped points. 

[0029] The graphical part and the internal part of the measurement scheme are bi-directionally linked such that the 
25 master-slave relation between the two parts is defined. The bi-directional link provides tfiat measurements can* be' 
' activated starting from the graphical part as well as starting from th o internal part. 
{0030]: This means that 

•• -either internal methods of the graphical objects maybe invoked by addressing the graphical content of the moas- 
30 . .urement scheme. After pointing to a geometric entity In the stencil, a method, appropriate for the type and location 
of the geometric entity, is invoked to map the entity in the image. For example, clicking on a line may invoke the 
user interface method to position two pairs of juxtaposed points, each point pair yielding a midpoint coincident on 
the requested line. In contrast to a passive bitmap display of the measurement objects, which would require a ' 
knowledgeable user to map them graphically, the measurement scheme supplies the appropriate method. Anal- 
35 ogously, a complete sub-graph corresponding to a measurement operator iii the rneasurement dependency graph 

may be invoked from the measurement template window. For example, when the mouse cursor is over a the 
graphical representation of a measurement operator in the template window, it may change from an arrowhead; 
appearance to a measurement operator symbol (e.g. a distance or an angle symbol), denoting that the associated 
sub-graph in the measurement dependency graph may be aetivatedby pressing the mouse cursor. 
to • or methods affecting the graphical behaviour of the measurement methods in the stencil and the.radiological image 
. may be steered by the internal methods or the intematcontrol flow* (e.g. a graphical entity in the stencf highlights 
when the corresponding node in the measurement dependency graph is visited). •••• . 

. [0031] In one embodiment a measurement scheme is selected from a measurement stencils repositajy. A measure-' 
45 - ment stencils repository comprises a number of measurement stencils. The selection is preferably performed iri cor- 
respondence with the examination type and/or the radiographic exposure conditions. 
£0032] The measurement scheme can be implemented as an active or as a passive, measurement scheme!. 
An active measurement scheme (also called measurement template) is differentiated from a passive one in that the. 
measurement procedure is.guided in part or completely by a computerised procedure operating on Image data. Guid- 
50.. ance may range from simple computer guided sequencing of measurement operations to fully automatic execution 
• based on automatically determining landmarks and deforming model contours representing anatomy. 
■ [0033] The graphical part of the measurement scheme can be implemented as an stencil-overlay on the' displayed 
X^y image. Thepolnts and lines of the scheme must then bedrag-ged by me user to me correct anatomical position 
in the displayed image. • '•';, 
ss [0034] fn one embodiment, display of the graphical part of the measurement scheme may be omitted when meas- : 
u'rements are simple and performed routinely on the radiograph. User guidance is effectuated here solely through the " 
measurement dependency graph, which may generate the measurement points and objects immediately in the image, 
after which the useris requested to map them to their actual position. An example of such simple measurement is the 
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cardio-thoracic index, which ratio is calculated routinely on thorax AP radiographs. 

[0035] In one embodiment, more than one radiographic image may be displayed and a measurement scheme may 
be activated with more than two associated graphical parts. This configuration is useful for performing 3D measure- 
ments from measurements on a limited number of projections. For example, 3D spinal measurements may be obtained 
from locations of anatomical points identified on frontal and lateral radiographs of the spine. A measurement template 
with two graphical parts is activated here; one for each of the projections, and each graphical part guides the user in 
the mapping of the projection point of a certain 3D point in the associated projection images. A measurement point 
node in the measurement dependency graph will compute its 3D coordinates from the locations of the point in the 
respective projections. Likewise, 3D measurement object and operator nodes have 3D methods to calculate the 3D 
parameters of the object and the 3D measurement values respectively from the values of their children nodes. 
[0036] Specific features and further embodiments of the present invention are set out in the dependent claims. 
[0037] One aspect of the present invention relates to a computer program product adapted to carry out the method 
of the present invention when run on a computer. The computer program product is commonly stored in a computer 
readable carrier medium such as a CD-ROM. Alternatively the computer program product takes the form of an electric ' 
signal and can be communicated to a user through electronic communication. 

[0038] The method of the present invention is advantageous over the prior art for the following reasons. 

A first advantage of using a measurement template is the unambiguous and instantaneous indication as to the 
position of the measurement objects in the actual image. That is, the first source of random efror in the prior art, 
due to a badly or non-defined measurand, is eradicated. 

In a film-based method or in a generic computerised measurement method in conjunction with a radiological meas- 
urement atlas, repeated focusing of attention back and forth between the atlas and the radiological image are 
needed to locate the points precisely. 

Moreover, in the absence of an atlas scheme, the position of measurement objects is not defined and hence 
different users may locate a given anatomical landmark differently. 

In accordance with the present invention, the use of a digital template, especially when displayed in proximity of 
the radiographic image, substantially reduces the effort to locate measurement object points precisely and hence 
reduces intra-observer measurement error (error in repeated measurements of the same quantity performed by 
a single user) and inter-observer measurement error (error in repeated measurements of the same quantity be 
tween different users). . 

In the present invention the measurements are digital in contrast to the analogue measurements of the filrrnbased 
method. The precision of positional location is thus dependent on the image resolution (physical size oft a pixel 
unit; which typically is around 100um for digital detectors). 

In a specific embodiment measurement errors may be further reduced by the possibility of zooming into portions 
of a radiograph prior to locating a measurement point, and occasionally applying a sub*pixel location algorithm to 
define the position of image points with precision ,up to a fraction of a pixel unit. : 

". 'Interesting measurement quantities are often den^fram basic measurements by calculations that combine, meaty 
measurements into final quantities, orfrom transfonriations, filters andfitprocedures applied on the measurements. 
In such cases the uncertainty characteristics of the derived quantities dan be derived by error propagation. 

: The ensembles of digital measurements laid out In a measurement s^emecoflectively define a standard manner 
,. on how the procedure must be performed. Standardisation is; advantageous because it enables objective compar- U 
ison (e.g. in second opinion gathering). , ; 

By the use of a measurement scheme the order of placement of measurement objects In the image is fixed and 
,. logical (eventually the ordermay be customised within the constraint or dependency). Therefore ( in the template- 
guided measurement method of the present invention,, thereis no need to study the measurement scheme in an 
;, atlas to.detect the order in which to place the measurement objectsso as to guarantee that all depending objects 
. t are placed^rst This way of operation is fyndamentelly Afferent frorn the film-based method, wfiere measurement . 
templates may even be absent altogether. Repeatability of measurement procedures, obtained by an imposed; 
order, enhances throughput of diagnostic evaluation of the radiograph, which is important for clinical departments 
such as emergency radiology and intensive care units! 



The automatic abnormal value signalling functionality provided by a specific embodiment of the method of this 
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invention, operating instantaneously on measurement values as they become available, -mikes-the necessity su- 
perfluous of haying the normative reference tables at hand, such as published in the open literature and meas- 
urement atlases. — 

• The method of the current invention is particularly suited to implement different measurement schemes in the field 
of paediatric radiology, since measurement templates can be customised to reflect changes in paediatric anatomy 
over time. 

The graphical part of a measurement template will display the measurement points, objects and entities using 
anatomical landmarks corresponding to the patient's age. Furthermore age-specific normative values associated 
with the measurement template can be applied. 

• In the field of orthopaedics the invention is particularly useful because instantiated measurement schemes can be 
stored and used for further follow-up after therapy. The graphical part of an older template can be retrieved frorp 
memory andean serve as a new template for a current measurement. The older measurements can be compared 
with the new measurement results to assess evolution of a patient's condition. No confusion as to the manner in 

- which the measurement objects and entities were defined is possible. Furthermore, second opinion gathering 
becomes more objective when it is based on the same graphically defined measurement scheme. 

• The method of the present invention is particularly useful in emergency radiology since the measurement scheme 
relies on the placement of a set of measurement point only after which the auto-completion of all dependent meas- 
urements is implied and automated so that the procedure lends itself to quick execution. 

This feature is not available in the prior art based on manual measurement based on analogue film measurements 
or computer measurements with a generic measurement tool. 

•• As it will be clearfrom the description in the sequel, both repeatability and reproducibility are substantially enhanced 
by the stencil-guided measurement method of the present invention. 

10039] Repeatability conditions include (a) the same measurement procedure, (b) the same observer, (c) the same 
measurement instrument, used under the same conditions, (d) the same location and (e) repetition over a short period 
oftime. 

Whereas fulfilment of these conditions is not guaranteed by the prior art methods, conditions (a), (c) and (d) are met 
by the current invention because use of a stencil and programmed methods to define the measurands ensure an 
identical procedure applied under all circumstances. Conditions (b) and (e) are fulfilled because the instantiated meas- 
urands in the image according to the stencil are stored: and supplied to other referring clinicians, who only need to 
perform a confirming analysis. 

J0040I A valid statement of reproducibility requires specification of the conditions changed. The changed conditions 
may include: (a) principle of measurement, (b) method of measurement, (c) observer, (d) measuring Instrument, (e) 
reference standard, (f) location, (g) conditions of use, (h) time. Whereas fulfilment is not guaranteed by prior art film- 
based methods and generic computerised measuring tools, the stencil-based method as laid out in tbecurrient invention 
enhances reproducibility considerably. In particular, reproducibility conditions (a), (b), (d), (e), (g> are met because the 
definition of the measurands and their determination are fixed by andiaid out In the stencil. The stencil-based method 
of the present invention is further invariant to reproducibility conditions (f) and (h), and condition (c) is achieved because 
storage and retrieval of prior instantiated measurement schemes render the need to re-perform the measurement 
scheme superfluous. 

(0041} The objective of a measurement fe td^^ 

quantity to be measured. A measurement therefore must begin with an appropriate specification of the measurand, 
the method of measurement, and the measurement procedure. In general, the result of the measurement is only an 
approximation or estimate of the value of the measurand and thus is complete only when accompanied by a statement 
of the uncertainty of that estimate. J - '•' 

Pridrart methodssuch as film-based methods and generic computerised measurements do not cope with the problem 
that neither measurands nor mieasuremerit method aw ejeflned.' 7 *'* '. ."" ••■'".'■"v: 

- As a result, the uncertainty of the estimated value of a .measurand cannot beidetermined because (a) a measurand 
must be defined by a standard method of measurement and (b) the implementation of the standard measurement 
method. . • • • •.•*•. 

in the steneiMguidedmeasurement method of the present invention, both definition of the measurands and implemen- 
tation of the measurement method are prescribed. 

£0042] Specific embodiments of the present invention wilf be explained with reference to the following drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS :^-T.-— ■ 

f0043] 

Fig. 1 generally illustrates the method of the present invention, 

Fig. 2 illustrates the manual mapping of measurement objects according to the geometric specification of the 
measurement scheme, 

Fig. 3 illustrates the ROI-magnified manual mapping of measurement objects according to the geometric specifi- 
- cation of the measurement scheme, 

Fig. 4 illustrates the automated placement based on deformable contour segmentation. 

DETAILED DESCRIPTION OF THE INVENTION 

[0044] Priorto the description of the flow of operations of the method of the present invention; four different imple- 
mentations of measurement scheme will be described: 

as a graphical object, 

as a tree or graph oriented graphical structure, 
- - as a-tree or graph oriented textual structure, 
as a structured document. 

• As a graphical object: 

£0045] In this representation the anatomy is schematically depicted as a collection of outlines of bone and other 
radiologicaliy well-manifested landmarks. The measurement objects are drawn superimposed onto the anatomical 
outlines. 

This representation is useful in the courseof performing complex measurement procedures since in this implementation 
the. position of the measurement points is displayed with respect to the anatomy of the patient. 
Fundamental to the representation as a graphical object is that this representation graphically depicts the;spatial re- 
lationship between anatomical features and measurement objects. • 

The main advantage of this view is that it enables the user to define user points unambiguously. Hence this represen- 
tation enhances the precision of the measurement procedure (i.e. the error that would result on the measurements 
when performed repeatedly by different radiographers). 

Significant spatial relationships for two dimensional (plane) images can be defined for example (analogous relations 
can hold in three or more dimensions) between points, a point and a line (more general a point arid a curve), a point 
and an area, between lines, line-area, between two areas, etc. Examples of relationships areincidence, betweenness, 
containment, etc. 

[0046] For example, the center line of the femur bone, which is a 3D volume, may be defined on the basis of two 
points, each of which are exactly in the middle between a pair of juxtaposed points on the femur contour. Thecircularly 
shaped femur head may be described by a circle segment fined only the femur head edge pixels (circle-contour inci- 
dence). 

' As a tree or graph oriented graphical structure . !; 

[0047J The main advantage of this representation is the ability to view the dependencies: between measurements, 
and thus the implicit order of performing the. actual measurements. 

Measurement entity sharing, the term used in this context to denote the use of a measurement entity by more than 
orie super entity,- is also expressed in the graphical staicture. 

When one or more entities are shared by more than one super entities, a dependency, graph is obtained, whereas in 
the absence of sharing, a dependency tree results . 

Obviously, measurements that depend on others can only beperformed or calculated when the latter have been done. 
•The procedural nature-implied by this representation is an.advantage because it guides the User and formalizes the 
measurement, . ■ . 

in. the prior art, it was the user's task to analyze ..and extract the right order from -a published measurement scheme, ■ 
wtiich is time-consuming and error-prone. In the absence of a guiding scheme, measurements performed at a given 
point In rime may become 'invalid When depending measurements were erroneously skipped over. 
[0048] A measurement scheme will typically have more measurement results as outcomes. Each measurement 
result is represented by a graph, or a tree when there are no cycles, i.e. no shared measurement entities. A collection 
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of disconnected trees is called a forest. - i— ™.--— 

[0049] Though the definition of a tree implies no direction on the edges, a tree is usually drawn with the root at the 
top and the direction is from top to. bottom. The root node represents the highest-level measurement entity, usually 
having a diagnostic meaning, a name according to the medical nomenclature and possibly an associated normative 
values table. The children of the root node represent constituent measurement entities such as points, lines, circles 
etc. Each of these nodes may in turn depend on still other measurement entities. The bottommost terminal points in a 
tree are called leaves. Because the most primitive user interaction by means of which anatomical points are indicated 
in the digital image is a mouse pointer click, all leave nodes graphically represent points (at pixel or sub-pixel locations) 
and all other graphical objects such as lines, circles or curves will ultimately be determined by user-defined points or 
calculated points. 

[0050] All nodes are collectively called measurement entities. Each node basically falls into one of two types: 
measurement operator 

The measurement operator performs a measurement operation. This type is further subdivided according to the 
type of resulting quantity. 

Pure measurement operations 

Pure measurement operations have measurement objects as their arguments and produce physical quantities 
such as distance, curvilinear length, angle, area, volume etc. as a result. 

For example, in radiography, the heart volume in a neonate is calculated as V ~ F ■ A B C, with A,B and C 
perpendicular axes of an ellipsoid inscribed in the cardiac silhouette on AP and profile radiographs, arid Fa 
factor including both the ellipsoid formula constant and the magnification factor, and further beingdependent 
on patient specific data such as age and weight. 
Arithmetic measurement operations 

Arithmetic expressions operate oh measurement objects and result in (clinical) indices such as ratios between 
distances (e.g. cardio-thoracic index), ratios between areas (e.g. ventricle-brain ratio) etc. 
For example, the cardio-thoracic ratio is defined as CTwith ethe width ofthe cardiac silhouette at its greatest 
extentto the right and left of the spine, and 7"the inner width of the thoracic rib cage measured at different levels, 
fn this category fall complex arithmetic relationships-expressed by a nomogram in the prior art. For example, 
the calculation of the body surface area and the assessment of gestational age from lumbar spine dimensions 
are based on nonlinear relationships of the constituent measurements or parameters. The graphical evaluation 
of the nomogram is emulated electronically in the present invention. 

measurement object 

A measurement object is a geometrical object that can be depicted graphically, and that has physbal or ana- 
tomical support fn informatics terms, measurement objects are instances of a class. 

The following classes are particularly relevant: . 

- Point 

A point represents a location in space, on the Image and on the stencil drawing, and has no width, height or 
depth. 

Measurements based on points must adopt a convention for locating the pixel's origin. Because each pixel 
has a certain size, the point represented by a pixel can be either the center of the pixel's square, or one of the 
comers of it 

The point class is further subdivided in 1 
User-defined points (Upoints) - 

These points must be drawn on the screen by the user either by a mouse pointer, or other digital input 
means such as a graphical tablet e.g. 
- Derived points (Dpoints) 

These points are needed In the course of the measurement procedure but in contrast to user-defined 
points, derived points are the : .result o,f geometrical calculation based on available graphical objects. Ex- 
amples are the intersection of two lines, the intersection of a line and the perpendicular of a line through 
a given point etc. The procedures to locate a certain point as a geometric locus are defined in the methods 
'of the point class. ^ . V ' ■ ' : 

■ ■ All' oth^ measurement objects in the digital image are ultimo 

derived. . • /*.••' 
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- Line - ^T r -r 

A line is a graphical object in 3D space with indefinite length and no breadth. A straight line is the shortest 
distance between two points in 3D space and is commonly referred to as "line". For twjj-dimensional images., 
we consider the projections of 3D lines in the image plane, which are also lines. In its simplest form, a line 
may be specified by drawing two points, each point drawn independently according to methods outlined in the 
point class. Alternatively, a line may be derived according to properties of other graphical objects (e.g. tangent 
and perpendicularorientation). Finally, a line may be calculated on the basis of a plurality of image points such 
as edge pixels by a fitting procedure. 

A line segment is specified by limiting the extent of a line between interval bounds. 

Compound linear objects such as triangles, quadrilaterals, regular and general polygons are classes based 
upon the line class. 

The line class specifies the methods how the line is determined in the image. 

- Analytic Curve ■' 

This class contains graphical objects that are generated in the image by evaluating an analytic equation. This 
class is also called parametric curve class . Two types of analytic curves deserve special attention in the conte>rt 
of a diagnostic measurement tool: 

• Circle 

A circle is a closed curve, all points of which are the same distance from a point called the center The 
center point may be drawn according to one of the methods specified in the point class. The radius may 
be specified by drawing another point on the circle's circumference. Also, three non-collinear points, each 
drawn independently according to methods outlined in the point class, uniquely define a circle. Alterna- 
tively, a circle may be specified according to properties of other graphical objects. Finally, analogous to a 
line, a circle may be calculated on the basis of a plurality (more than 3) of image points, such as edge 
pixels, by a fitting procedure. 

A circular areis specified by limiting. the circumference of the circle between angular bounds. 

♦ Ellipse 

An ellipse is generated by a point moving so that the sum of its distances from two points (the foci} is ' 
- - constant and equal to the major axis. The ellipse object is important in planar images because the parallel 

projection of a circle in 3D space is an ellipse. Metallic rings (e.g. embedded in the liner of a sphere- 
shaped cup prosthesis to make its outline and orientation visible in the image) with acertain 3D orientation 
project in the 2D image plane as ellipses, 

Again, as for all graphical objects having analytic description, ellipses may be calculated on the basis of 
a plurality of image points such as edge pixels by a fitting procedure. 

* Polynomial curves 

Polynomial: curves such as a parabola have clinical utility in measurement schemes as well. Forexample, 
the outline of corresponding; anatomical landmarks on the spinal process may be fitted with a polynomial 
of degree sufficient to make a least-squares fit error lower than a pre-defined threshold. Pathological 
displacements with respect to the mean curve of one or more vertebrae may then be detected by imposing 
a suitable threshold on the standard deviation (typically 2...3 SD). 

Pseudospondylolisthesi's and spondylolisthesis, which are manifested by a. forward shift or bulge of a 
vertebra, over the edge of the one below it, may be. assessed likewise. A curvature measure, computed 
on, the basis ot a fitted polynomial; may foe used to quantifykyphosis or lorclosis spinal deformity. 

• Bezfcr and Splrnc Curve -.- 

Curves in 2D and, surfaces in 3D with a more complex shape rfiay foe descrlbedby plecewise polynomial 
.: .. subdivisions of a desired cbject shape, each of the curve segmen ts defined by a control polygon and a 
basis of interpolating. Wending functions to establish the relationship between the eutve and the control 
. .polygon. TiiecontiT>rpolyo>nfo^ 

, c ^ .by met^ds outltned above.or may bederived on the basis of least-squares fitting. Of particular interest . 
v. : ; . ,is the approximation of an object shape in the image with a fitting BeztSr or spline shape based on the '' 
; = -edge points , of the object. Distance of an open-ended curve may be computed ,as the curvilinear path - 

length defined a„ lite path length mtcgt v between 2 points tying on the curve. .An angle between two lines 
tangent to the .curve at a two- points , may likewise foe calculated based oh the derivatives' of the curve's 
. analytic function. . ■ : : 

[0051] Because the result of a measurement scheme is a quantity, the Toot node type is always a measurement 
operator (either pure or arithmetic). • .... . - 
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[0052] The arity of the nodes of a tree will generally depend on the complexity of the measurement-Measurement 
operators such as distance or angle involve two operands (e.g. the distance between two points, between a point and 
a line, the angle between two lines etc.). Arithmetic expressions based on measurements can. also be expressed by 
or converted to a tree having binary nodes only. These binary arithmetic trees are also called aparse trees, and are 
constructed from expressions in either parenthesized infix or postfix (reverse Polish) notation. 
[0053] The arcs between a node and its descendants (children) denote the relationship between them. For meas- 
urement operators, the relationship refines the arguments; for measurement objects, relationship can be read from the 
descendant node (children) towards the depending node (parent). For example: 

■ the binary measurement operator "angle" subtended between two lines may be described as: 
angle (from line, to line) and the relationships are from (first line), to (second line) 

■ the intersection point between two lines is described as the point that is the geometric containment of both lines: 
intersection (containment line, containment line) ... f 
The term "intersection" is used here instead of the more general object name "point". The relationship can be read! 
the line contains the point. 

■ a circle may be defined by three point coincident on it: 

circle (incident point, incident point, incident point). The relationship can be read: the' point is incident on the circle. 

■ the perpendicular distance from a given point to a given line is computed by first calculating the intersection between 
the given line and the perpendicular onto that line, and next calculating the distance between the given poinl and 
the intersection point 

■ ratio (nominatormeasurement operator, denominatormeasurement operator) and the relationships are nominator, 
denominator: ■ . - . 

ratio (nominator operator, denominator operator), and each of the arguments is the result of a measurement op- 
eration ■ ■ ■ . 

[0054] Specification of the method of computing the object of a node may be intrinsically or extrinsically. The leave 
nodes, representing points in the image, are generally computed by a method based on internal parameters of the 
nodes only. In passive measurement templates, for example, the points may be defined by pixels designated by the 
user via the graphical user interface. 

Definition of the relationship enables a node to perform its computation of the value of the node. In the example given, 
when two lines each have a relationship of containment of a point, the method of computation of the intersection point 
is invoked and the result is passed to the higher-level nodes. When the relationships between a point and a line with 
a desired line are one of incidence and perpendicularity respectively, an algorithm that computes of the perpendicular 
line through a.point onto a line is performed. 
. For active measurement templates, the leave nodes of the measurement dependency graph represent groupings of 
pixels, belonging to an anatomical contour, and the method of computing these pixels is specified by an edge detection 
and linking algorithm. The depending node will subsequently perform a model-based deformable contoursegmentation 
oh the basis of edge and region based features, mapping the instantiated model contour onto the anatomical data. 
The output of such a node is one or more points, representing requested measurement point(s) mapped to its corre- 
sponding position in the image. 

* As a tree or graph oriented textual structure 

[0055] This representation is particularly Suited to perform customfzation of existing measurements, or to create 
measurement schemes from scratch. A two-pane view can be used for this purpose, ftie left pane displays the trees 
" making up the measurement scheme(s). To convey dependency structure, each sub-level of the tree is Indented with 
respect to parent level. Each node of the tree- has an associated flag indicating that the underlying sub-tree is either 
contracted (flag is a V sign, indicating that the node has children that currently- are invisible but may be viewed by 
-pressing the flag button, upon which the flag turns into a "-" sign) - or expanded (flag Is a "-" sign, indicating that the 
. node's children may be hidden by pressing the flag button, after which the flag turns into a V sign). Customization 
proceeds by copying sul vets of measurements from an existing scheme (or a numberof schemes) to a new scheme 
by mouse-dragging the top^-node of a sub-tree from the left pane to the right pane. Alternatively, a copy-paste sequence 
(CTRL-C/ CTRL -V key strokes) may be used. To refine the resulting measurement tree(s) in the right pane, individual 
measurements may either be deleted, or added by drawing them from a generic measurement toolbox. 

" As a structured document * • • . ' ■ 



-fOOSej lrt this representation, individual measurements, collections of measurements-andmeasurernenf dependen- 
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cies are stored in a document comprised of information entities. Layering type and structure^Bie-iftformation is en- 
coded according to a generally adopted standard such as the XML (Extensible Mark-up Language) recommendation 
for example. The main advantage of this type of information storage is the universal exchangeability of measurement 
schemes and associated measurement data. The declarative nature implied by this representation is advantageous 
for subsequent operations such as storage and archiving, abnormal value signaling, and data mining. The structured 
document representation also lends itself for customization of templates. Custom-made measurement schemes may 
be built from scratch, by performing measurement drawn from the general measurements toolbox, and finally saving 
the resulting procedure. Likewise, an existing measurement stencil can be modified by adding measurements from 
the general measurements toolbox or omitting measurements (both subject to conditional dependencies) and saving 
the result as a new template. A third way of customization consists in that a new template may be constructed departing 
from two existing complete measurement templates by matching and merging their XML trees (or sub-trees) according 
to a set of merging rules operating on individual nodes such as updating, adding, deleting and moving nodes or ag- 
gregated nodes (sub-trees) such as sub-tree graft (inserting a sub-tree), sub-tree prune (deleting a sub-tree), sub-tree ( 
move (moving a sub-tree). Node equality comparison may be based on matching of node ID's or the node tag name.' 
and content. Prior art for merging trees may be used to this purpose such as the IBM Alphawbrks' XMLTreeDiff and 
XMLDiff ("XML Diff and Merge Tool"). More than two templates may be merged by iterating this procedure. 
[0057] The following is a description of the flow of operations when passive measurement templates are used. 
[0058] After the digital image has been acquired and displayed and the measurement template has been selected, 
a geometrical mapping process is invoked to request the user to place all measurement objects in the image whereupon 
the corresponding position in the image is entered into the computing system. 
^[0059] Depending on whether or not image data or image-derived data are used in the process of locating the meas- 
urement objects, a measurement template .is termed either active or passive. 

[0060} In a passive measurement template all measurement objects are placed manually in the displayed image. It 
is the user's responsibility to locate all measurement points in accordance with the measurement template sketch 
displayed simultaneously and in close vicinity with (or superimposed on) the radiograph. 

[0061] When using passive measurement templates, two embodiments of user-system interaction may be consid- 
ered depending on the way in which the internal informatics model of the measurement template is activated. ~.-<v. 
[0062] In one embodiment, all user-requested measurement points are mapped prior to generating the measurement 
objects that depend on them. All measurement objects are highlighted on the digital sketch upon which a copy of the 
object is generated and the user is requested to drag and adjust the copy to its corresponding position in the actual 
radiographic image. When all measurement points are available, the depending measurement objects are generated 
and the result of the measurement operators is computed. -.- 
[0063] This mode has the advantage that it enables quick execution of the measurement template because all user- 
interactions are concentrated in a single pass. 

It has the disadvantage that fewer insight is gained how the template has built up the measurements. Afso, eventual 
«rore in locating Individual points are less easily discerned and •corrected, 

[0064] In a secondembodiment, the measurement objects are generated when all points defining them are available, 

and likewise, the result of measurement operators may be computed when all required objects are known. 

This mode has the. advantage that the userhas close control oyer the. measurement procedure as it is constructed 

gradually, which is especially advantageous for complex measurement schemes, In.the occurring event that the user . 

■establishes the current object wrongly; he/she may re-position the underlying points, without affecting the rest of the • 

procedure. 

It has the drawback that superposition of too many lines and mark-up in the course of the proceduremay prevent the .. . 
user from i accurately positioning the ; remaining points. . , . 
{006SJ W^enameasur^ento^er^ 

me measurement value is calculated. The calculated date may be corrected qn the basis of calibration data derived 
from calibration objects in the imago (such as a ruler). ■■ =. - - : •-. \ 

(0066} Furthermore, based on the measurement dependency graph, errors can be propagated to provide an estimate ; 
of measurement uncertainty. 

■f0067] . At the time of processing ameasurement operator node, a process can be iniflatedthat selects the normative . 
value associated with the anatomical quantity produced.by the measurement operator node in the given template: 
The proper identifier and patient data are sent to a normative values repository, and the raw normative reference and: 
its semantics are retrieved, 

Nomialisalion is performed when needed r for exu;r. r \, to express the acceptable range fo the commonly used Stigma 
range, or to a user customised range When the standard deviation is retrieved, or to perform the nomographic look-up 
procedure when normative values are encoded in nomograms. The final normative Values are applied to the result of '. 
the measurement operator, and nqrmafity/a^nomiality is messaged. 

[00681 The result of measurement operators and their associated uncertainty bounds can bevdlsplayed toxtually in 
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a measurement values window, and depicted graphically using suitable mark-up produced by anaatomatie dimension- 
ing process. Likewise, normal values can also be displayed textually and depicted graphically, to enable clinical as- 
sessment of the normality/abnormality of the measurement value in relationship to its uncertainty. 
[0069] Two modes may be considered when using a passive measurement template: manual placement mode and 
enhanced placement mode. 

[0070] In the manual placement mode all constituent points required to derive a given measurement object are placed 
manually in the image by dragging the mouse pointer until it is over the anatomical location of interest and clicking the 
mouse or a key to fix the position. Portions of the image may be zoomed out to define the point's position to within 
pixel unit precision. 

[0071] The enhanced placement mode is achieved by hinting the user as to the position of a constituent point of a 
measurement object by confining the placement to a set of points or objects defined by the graphical construction 
defined so far. These set of geometric objects are called the geometric loci of a geometric problem. 
. When the mouse cursor is moved, that point of said set, which is closest to the mouse cursor, is highlighted along with 
its property with respect to the measurement object being drawn (e.g. it is a tangent point, perpendicular point, center 
point, end point, topmost point, leftmost point, inflection point etc.). 

If the user accepts the suggested point, he of she presses the mouse button, even when the cursor is not exactly over 
that point, and the mouse cursor will immediately jump to it. 

This mode is also called snap placement mode, and makes it possible to exactly pick out the location of a given point. 
Obviously, because analytical expressions of measurement objects enable exact geometrical solutions, the geometrical 
measurement precision is substantially enhanced in this mode. 

Setting snap to lines or other graphical entities has a similar effect: e.g. the line tangent to a circle (there are two such 
lines) through a given point is selected when moving the cursor nearest towards the intended tangent point, and a 
mouse click will teleportthecursorontothattangentpointon the circle, after which drawing of the tangent line completes 
the drawing of the line measurement object. 

When the set of attractor points contains only one point, the drawing action completes without any further user inter- 
action. 

This auto^completion property is made possible because the measurement template names each measurement object 
(according to medical nomenclature) , and the measurement scheme uses the name of objects to uniquely designate 
each of them. 

-[0072] For example, suppose that three lines have been drawn in the image, named A, Sand C, and a fourth line D 
needs be drawn through a point p perpendicular to line B: 

Line Dean be drawn immediately after placing point p, because of all three possible perpendicular intersection points, 
the Intersection point with line B is the target point. 

Without naming of objects, the user would be faced with three possibilities in this example. 

When multiple snap objects still exist, such as in the case of a line through a given point and tangent to a circie, the 
user must resolve the intended point by moving the cursor towards it (which will invoke highlighting of it) and completing; 
the editing by pressing the mouse. 

Enhanced measurement precision and faster measurement completion are the main advantages of this mode 
[0073] Defining objects based on attractor or snap objects is different from active measurement objects, to be defined 
further, because they are geometric loci, based on andcaJctifatedwith the interna! descriptions of the graphical meas- ; 
urement objects (e.g. equations of line, circle...) whereas active measurement objects get tneirfinal position and shape: 
based on Image or image-derived data. 

[0074]. The following is a description of the flow of operations when active measurement templates are used. - 
[0075] Also inthls embodimentthe digital image that has been acquired lis displayed. Then a measurement template 
is selected anda geometrical mapping process Is Invoked to request the user to place admeasurement 
Image whereupon the corresponding position In the Image is entered into the computing system. 
When image data and/or image-derived data are used in the process of locating the measurement objects, the meas- 
urement template is termed active, 

[0076] In contrast to the passive measurement template, the geometrical mapping process for active measurement 
templates is based on automatic landmark, feature and object locating by applying suitable image processing operators 
and algorithms. '■ • ' -' ■ 

[0077] These algorithms assume that the location of the -measurement points may be modelled either directly or 
indirectly. ' v^;;;Y- ;: : : 

In direct/modelling the measurement point is specified in terms of a local pre-defined configuration of grey value char- 
acteristics/in indirect modelling it is assumed that a set of measurement points may be defined in spatial relationship 
: to an anatomical shape of which representative models (either 2D or 30) have been built that are mapped in the image" ' 
(either 20 or 3D). In 2D Images, models may fypicafly take the form of collections of open-ended or closed curves.' In ; 
3D Images the models may take the form of surfaces: • : ': 
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[0078] Initialisation comprises the mapping of anchor points in spatial relation to which"tbieTrHi»» position of the 
instantiated model contours will be positioned. 

The initial pose and shape of the model contours is adapted according to the spatial disposition of the anchor points 
in the image. 

Pose comprises the translation, orientation and size (i.e. position of the origin point of the instantiated model contour, 
orientation of an axis of the contour relative to the orientation of lines connecting the anchor points, and the applied 
magnification). 

Shape comprises size and initial form of the contour, which may be expressed by a transformation matrix calculated 
on the basis of the anchor points of the model space and the corresponding anchor points located in the image. 
[0079] For example, two points in a plane define an Euclidean similarity transformation (translation, rotation, scale) 
in that plane, three anchor points (forming a triangle) define an affine frame (translation, rotation, scale and x and y 
shear), four anchor points (forming a quadrangle) define a projective frame (translation, rotation, scale, x and y shear 
and projections). 

On the basis of the computed transformation parameters, the initial curves representing anatomical entities are instan- ' 
tiated in the image, or the regions of interest around anatomical landmarks are mapped in the image. 
Initialisation further comprises setting suitable parameters for the segmentation such as parameters for the internal 
and external forces of snake deformation (i.e the coefficients of elasticity, rigidity and the weight associated with the 
image data). 

[0080] An image processing procedure subsequently calculates the precise segmentation of anatomical entity con- 
toure.or anatomical landmarks on the basis of image data only 

[0081] Measurement objects, the definition of which is determined in relationship to the delineated anatomical ge- 
ometry, are computed and drawn superimposed on the segmentation. 

The user is prompted to accept the segmentation and landmarks, and the computed measurement objects. In the event . 
of mal-positioned objects, the procedure is re-iterated with user-adapted initial parameters. 

[0082] The remainder of the flow of operations in an active measurement template is identical to the flow based" on 
a passive measurement template. 

[0083] Mixed measurement schemes may be constructed and effectuated also, where measurement, objects sare; 
located by combining manual mapping of the anatomical points with automatic point mapping based on contour, and 
region segmentation of some anatomical entities that have linked the required points into their representation. 
[0084] The following is a general overall description of an embodiment of a procedure for performing a digital stencil 
guided measurement according to the present invention. The constituent method steps are explained in detail. 

• Image acquisition and display 

[0085] A digital image representation of an image onto which geometric measurements are to be made is acquired. 
[0086] A large variety of image acquisition systems are applicable. The radiation image can for example be recorded 
on a photo-stimulable phosphor screen. The screen carrying the radiation image is then read put by scanning, it with 
stimulating radiation, detecting the image-modulated, light which is. emitted upon stimulation and by converting. the, 
imagerwise modulated light into a signal representation of the radiation image. 

[00871 In an alternative embodiment the radiation image is stored in a direct radiation sensor rendering [ a signal: . 

[0088] Stiff further alternatives for acquiring a digital signal representation of a radiation image, such as scanning of 
radiographic film carrying a radiation image, are possible. 

[0089] It wilt be clear that the image representation may also be retrieved from art imagedatabasej e^- part of a . 
RIS (radiology information system) or HIS (hospital information system) or PACS (picture archiving and communication 
system). .■ ■ ." . . . V 

[0090J The acquired image representation is next applied to a display device cpnnectedto acomputerand the image 
is displayed. 

♦Template retrieval and display. , . ... . . 

[00911 The measurement template is selected from a measurement stencil repository stored in the computer, 
[0092] . The selection is preferably made according to the examination type, imaging parameters and patienUpocifio . 
information. "* '.' ' .. ."" - ■ ■ • ■ , „ • 

The examination type determines the anatomical part imaged. Imagingparameters such as angle of jnctdence,.distance ; ... 
between radiation sou rce and patient and resolution are important selection.crfterions. Age and gender may afeQaffeet * 
anatomical appearance of body organs. ... .:• 

[0093] The selected template is retrieved, loaded and activated in the system. 
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[0094] In one embodiment this involves: 

- Activation of the measurement dependency graph and registration of all needed measucement methods for the 
constituent measurement points and measurement objects, 

- Generating a graphical representation of the measurement template displayed in close vicinity of the radiographic 
image. Associated with the graphical representation, a structured document frame is generated. The document 
frame holds a description of the graphical layout of the measurement scheme in the current image in a universally 
exchangeable format such as SVG (Scalable Vector Graphics, which is a language for describing two-dimensional 
vector and mixed vector/raster graphics in XML). In this document frame all actualised graphical parameters of 
the measurement objects and entities will be stored. 

Generating the measurement values window displaying the medical nomenclature of all measurands in the 
scheme, and a set of numerical values pertaining to the measurement. This set comprises the actual measurement 
value and units, an applicable calibration factor, the normative reference value if one is provided, and an abnor- 
mality-alerting box. Associated with the measurement values window, a structured document frame is generated 
holding a tagged description of the measurement scheme in a universally exchangeable format such as XML 
(Extensible Mark-up Language), and into which all measurand-specific measurement values will be inserted. 

• Locating measurement points 

[0095] The measurement template has the purpose of defining the measurands with respect to the anatomy and of 
defining how the constituent points, defining a measurand, are to be placed in relationship to the anatomy. 
[0096] In the described embodiment of the present Invention, three placement modes are described: manual place- 
ment, ROI magnified placement, automated placement based on deformable contour segmentation. 

1 . Manual placement (fig. 2) ' " 

[0097] In this mode, alf points are placed manually in the image by cursor clicks. 

Irvthls embodiment the measurement template imposes the placement order by highlighting each point.in sequence 
(e.g. by blinking the point in the template). 

Following this order ensures-that at any point in time, the system knows which measurement object is being defined 
by virtue of its definition in the template, and hence which measurement objects and mark-up are to be generated next 
or which measurement entities may be computed next because of availability of all constituent points of the measurand. 
This auto-completion feature is possible because the dependency graph, of the measurement template defines how 
measurement objects and entities are computed on the basis of existing objects. 
This property is not available in generic measurement tools of the prior art. 

[0098J Alternatively, completion of objects and entities of higher layers may be deferred until the user determined alf 
points of the bottommost layer. This alternative embodiment has the advantage that no graphical display of auto- 
generated objects and entities will obscure the precise placement of subsequent points, which is particularly useful 
when complex measurement schemes need be fjerfonmed. 

{0099J Alternatively, image data and measurement mark-up may be organised and stored jn different layers Into, 
computer memory, the display of which is turned on and off depending on the current action being performed. For 
example, when the user is positioning points, the display of anjr unimportant graphics may be turned off; 

2. FtOI-magnified placement (fig. 3) 

[01001 J" mode, measurement points are placed more accurately In the Image by magnifying a region of interest 
(ROI) box around the measurement point. 

Each measurement point has an associated ROI, which is an image region covering a particular body part or joint, and 
wherein the measurement point is contained. 

A ROI may contain a set of measurement points If all of them have similar anatomical semantics and are optimally 
detetmined simultaneously. 

For example, one single ROI covering the femur head is tised to locate three separate points sampled identically and 
approximately equidistant aiongthe femur head's contour and collectively approximating '(he circularly shaped femur 
head boundary. ; /■ * 

[0101] In the described embodiment, in order to determine the image ROI corresponding to a moasurcmen; point, 
the portion of the Image around the mouse cursor is magnified constantly as the mduse is moved over the image, the' 
size of said image ROI being determined by the measurement template. A mouse cursor dtekfreezes ; the magnification 
around the current mouse position when me operator believes the ROf is covering the image region containing the 
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desired measurement point. The point editing is completed by a second mouse click at the anatomies! image location 
of the measurement point As the magnification enables to view image data to pixel-level resolution, individual pixels 
may be selected unambiguously. _ 
[0102] An alternative for speeding up the determination of the image ROI's is the semi -automated position of ROI's 
relative to a few user-defined points, the position of which is anchored with respect to the anatomy as depicted in the 
measurement template. 

[0103] Anchor points are sets of points in the template, which have specific anatomical meaning, which are easily 
locatable, and which are few in number yet enable full location of all ROI's. 

In a full-leg examination, eligible points are the bony landmarks corresponding to the greater trochanter, the femoral 
intercondylar notch and the upper ankle joint center. 

Because each ROI has sufficient extent and its location around the measurement point in the image allows certain 
positional variation, mild positional variation of the anchor points is tolerated as well. The position of previous meas- 
urement points inside their associated ROI's may refine the position of ROI's of subsequent measurement points-. . 
[0104] A portion of the graphical measurementstencil around the measurement point maybe magnified also to better' 
view the point's relationship to the depicted anatomy. The measurement template encodes the precise boundary of all- 
ROI's around their measurement point(s), 

[0105] The zoomed portion of the image corresponding to the ROI may be placed either inside the image, at a position 
centered around the mouse cursor, or alternatively in a separate zoomed image data window. 
[0106] Analogous to the manual placement mode, generating the graphics and mark-up of measurement objects 
and entities may be deferred until the position of all points of the bottommost layer of the dependency graph of the' 
measurement template have been determined in the image. 

3. Automated placement based on deformable contour segmentation (fig. 4) 

[0107] In this mode, no direct placement of measurement points is required. Instead, the position of measurement 
points is defined in relationship to contour models representing patient anatomical entities. 

[01 08] Hence, the position of measurement points and objects Is Indirectly specified by .their associated position and 
local features on an instantiated and subsequently deformed model contour. 

Examples of the relationship of desired measurement points may be incidence (of the measurement point ortthe 
- contour), midpoint (with respect to pairs of specific locations on the contour), and barycenter (with respect to configu- 
rations of specific point locations on the contour). 

Examples of the relationship of desired measurement lines may be fitted tangent line (with respect to a specific location 
on the contour), symmetry line (with respect to juxtaposed segments of contour pixels). 

Examples of the relationship of desired measurement curves may be fitted circle arcs (with respectto specific circularly 
shaped segments of the contour), and fitted parabola's or higher order polynomials to extremal points of. groups of 
contours. 

Automatic placement and deformation of the instantiated model contour(s) can be based on automatically, segmenting 
the bony outlines based on Point Distribution Models (PDM) and Active Shape Models (ASM). This segmentation 
principle can for example be applied to individual .vertebrae of the spine,, the hand and the hip. 
[0103) . The initial position of the instantiated model contours can likewise be based on a small set . of anchor points,: . 
relative to which the initial pose (translation, rotation, scale) of. each instantiated contour is determined and further be 
adjusted according to affine or projective transforms. Alternatively, object recognition techniques maybe.used to locate- 
the anatomical site corresponding tq the model contour. 

* Chaining of measurement operations; . ~, . . <. 

[01 10] The computation of measurement objects and measurement entities follows an execution How as determined 
by the order of visiting the nodes of the dependency graph. ' . 

The bottommost level contains points only, of wbich tJie determination of their.position is requested from.the user. 
Levels above the bottommost layer contain measurement objects and entities that build upon these user points. 
[0111 J Prior art is used to traverse the measurement dependency graph. Because all measurement points must be 
determined first, prior to being able to execute the parent nodes, the visiting order is: . . . 

"visit the left subtree, then visit the right subtree, then visit the root" , or equivalent^, "visit the children, then their parent 
<s)«. 

This is a postorder graph traversal, based on a stack data structure to store the nodes. From a point of view of calcu- 
lation, this is a demand-driven data flow network, because a node can only be computed when all data, required in its - 
computation, are available. 
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• Derivation and drawing of measurands and measurement mark-up . - ^ - 

[01 12] Based on the measurement point's coordinates in the image, the measurement objects (lines, circles, parab- 
ola's etc.) are generated using methods that are widely known in computer graphics. 

• Calculation and display of measurement values and uncertainties 

[01 13] The calculation of the measurement operator values, such as distances and angles, are based on the geo- 
metrical location of the measurement objects in the image and are computed following principles of geometry. 
[0114] In this embodiment the raw value obtained from the image-based location of the measurement objects is 
corrected using calibration factors. 

E.g. magnification and pixel size are basic calibration factors to convert pixel distances to physical distance. 
Raw values, correction and calibration factors, and physical measurement values are tabulated and displayed in a 
measurement values window. ;' 
[01 15] Measurement entities such as linear and angulardimensions are represented in the radiographic image using 
graphical mark-up. Dimensions include dimension text, dimension lines, extension lines, and arrowheads. Type and 
appearance of dimensions may be customised to user preferences. 

[01 16] Measurement values may be displayed either discretely or continuously in the measurement values window 
as the user moves the position of a point over the image. 

In the discrete mode, the values are updated after user control is released (e.g. the user stops dragging the mouse 
cursor and releases the mouse button when he/she believes the measurement point is well located in the image). 
In the continuous mode, the values are adapted continuously during moving the measurement point in the image. This 
mode has the advantage that the abnormality alerting mechanism, described in the sequel, is continuously calculating 
whether some measurands are outside their normal range in response to repeatedly moving the measurement point 
around. 

When a measurand is on the threshold of abnormality, the abnormality signal will start toggling, alerting the physician 
that care must be exercised, and more confidence will be gained in the outcome of the measurement when finally the 
mouse cursor control is released on a considered position in the image. 

[01 1 7] Measurement results have some inherent uncertainty because the measured result is only an estimate of the 
■ "true" or'correct" value. To establish this relationship between the measurement and the true value of the measurand, 
the contributors to uncertainty in the measurement are established. In most cases a measurand Vis not measured 
directly, but is determined from a number N of other quantities X^,X 2 ...X N through a functional relationship f : Y=i{X A , 
X 2 ...X N ). The input quantities Xi,Xz,..X N upon which the output quantity ^depends may themselves be viewed as 
measurands and may in turn depend on other quantities as expressed by the measurement dependency graph. In the 
context of the present invention, the measurands and depending measurands are based on the determination of points 
in the image, and therefore, the measurement uncertainty will be a complex relationship /of the uncertainty in pixel 
location. An estimate of the measurand Y, denoted by y, is obtained from the functional relationship fusing input 
estimates x v x 2 —x N for the values of the rV quantities Jf, ,X 2 .. Thus the output estimate y, which is the result of the 
measurement, is given by y = fa ,* 2 ...X/v). Each input estimate x> and its associated standard uncertainty u(xj) are 
obtained from a distribution of possible values of the input quantity X^ This probability distribution may be either fre- 
quency based, that is, obtained from a series of k observations of or it may an a priory distribution. These types* 
are referred to as Type A or Type B evaluations of standard uncertainty respectively. For example, assuming a linear 
resolution 8x of the image, the value of a pixel determination process that produces a given location X"can tie with 
equal probability anywhere in the interval X - 5x to X + 8x and the process is described by a rectangular probability 
distribution of width Sx with variance ifl =(6x) 2 /12lmplying a standard uncertainty of u=0i29«x forthe linear position of 
any plxef in the Image (referred to as Type B evaluation of standard uncertainty). The standard uncertainly of y, where 
y is the estimate of the measurand Kand. thus the result of the measurement, fe obtained by appropriately combining 
the standard uncertainties of the input estimates x^..*^ This precis is cafed error propagation. 
[01 181 The square of the combined uncertainty u^y) for independently assumed input quantities and sufficient lin- 
earity of fat the Xf is given by 




me partial derivatives to X, being evaluated at ^called the sensitivity coefficients. For correlated input quantities, the 
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expression for the combined variance i?(y) associated with the result of a measurement is- -- 



M N N ^, 3 , 



with i^XfrXfi^uixfruix^XrXj) the estimated covariance associated with x,and X/ and the degree of correlation between 
xyand Xycharacterised by the estimated correlation coefficient rix^xj). Although uj,y) can universally be used to express 
the uncertainty of a measurement result, in the context of a medical application, a measure of uncertainty is given that 
may be expected to encompass a large fraction of the distribution of values that could reasonably be attributed to the 
measurand. The expanded uncertainty denoted by U is computed by multiplying the combined standard uncertainty, 
u c (y)by a coverage factor k : U = ku c (y) and the result of the measurement is expressed in the measurement values ' 
window as Y=y±U. Usually, k is set to 2, meaning the uncertainty is denoted at a 95% (2SD) confidence level. ; 
[0119] The uncertainty bounds of a measurand may also be depicted graphically in the image by additional lines 
drawn at either side of the measurement mark-up line of the measured quantity in its associated graphical, layer. 

• Retrieval of normal values and comparison with measurement values 

[0120] After selectiofv-of a measurement scheme from the measurement stencil repository, the associated normal 
values table for the measurement scheme is retrieved from the normal values repository. 

When a measurement entity node of the measurement scheme is visited in the measurement dependency -graph, the 
medical nomenclature of the measurement entity and patient-specific data (such as gender, age, length.:.} are used 
as a selector in the normative values table of the measurement scheme. ■ . 

The normative value, and the associated comparison function are retrieved from the normative values table, after which 
measurement value and normative value are compared. ■„..-' 
An abnormality signalling function alerts the physician should the measurement value lie outside the normal range. 
This functionality is instantaneous, that is, normality or abnormality of a measurand is signalled immediately when the 
result of the measurement value computation is known. This ability is impossible in the prior art, where themeasurement 
value needs be looked up afterwards in the proper atlas,,. 

[0121 J Normal values and other' numeric information useful for measurements in radiology are calculated and-ipre- 
sented in many different ways, for example as 

- a fixed cut-off value. 
Examples: 

- any value of the angle of the epiphyseal axis with the cylindrical bone axis different from 0° may be indicative 
of a metaphyseal fracture .. 
. - a metacarpal If index (which is 1 00 times the sum of cortical thickness divided by width) smaller than 44 may 
be Jndfeatn/e of osteoporosis 

- • a demographic and patient data-specific cut-off value. This case differs from the previous oae.fn that thocut-off • 

: pr threshold value- is made dopendem-on e.g. the age and gender of the patient, efhie race, body weight, height, = ; 
body surface etc. Hence, the cut-off value Is a function of one or more variables, which may bo stored in discrete 
form by sampling the function at suitable Intervals. The electronic calculation re-samples -fne interpolated function.-. 
. to retrieve thefjnal threshold^ ■ .. .. , 

- a normative Value presentedwith statistical, upper and lower limits, these limits.being defined as the mean value 

- plus or minus a number of standard deviations (SD). This latter tolerance: range, is normally equal to -2.SD to +2SD, . . 
. and means that about 2% of a normal population Willbe assessed as abnormally large and around 2% as abnor- '. 
malfy small with respect to the parameter assessed. As with the foregoing case, the tolerance range may be 
function of a number of variables. 
■ a threshold based on a multiple of the standard deviation, the SD being calculated on the basis of the deviation . : 
with respect to an ideal fitted geometric object (fit line, fit parabola, fit circle, fitted ellipse etc.) Example: vertebral,- - 
displacements due to dislocations and diagnosed by. means of.the perpendicular distance of anatomical landmarks 
with respect to an analytic curve fitted through, corresponding anatomical.landmariaaf a:sequenee of vertebrae. " . 
a nomogram (or nomograph), which is a graph containing several (usually three) parallel. scales graduated fof. . 
. different variables so that when a straight line connects values of any two, the related value may be read directly 
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from the third at a point intersected by the line. It assists in estimating data that normally wouttr require intricate 
calculations. For example, carpal length relationships are usually determined by a nomogram, where the grade of 
deviance from the mean for a specif ic child is established by placing a ruler between the two points on the scales 
that correspond to specific measures in the hand radiograph of the children. The electronic implementation stores 
the axes in discrete form. To retrieve the grade of deviance, the position along the axis is interpolated from the 
value of the known input variables, followed by an emulation of the graphical construction to retrieve the requested 
deviance given the input variables. 

[0122] The measurement template provides a coupling between an individual entity and its normative value through 
a link to a record in a normative reference database. 

The normative reference database is organised around tables, grouping mean value and standard deviations for each 
measurement scheme in one table. 

Tables may further be organised hierarchically according to a whole-part relationship (e.g. a normative table for spinaj 
measurements maybe composed of tables forthe individual vertebrae). The link may be through the generally accepted 
medical name for the anatomical quantity or it may be based on a unique identifier when the quantity is referred to by 
a plurality of names according to different translations or user preferences. 

[01 23] Just as measurement values and their uncertainty bounds are displayed graphically, the normal values may 
be depicted by suitable graphical mark-up as well. For example, along with the actual measurement value, the mean 
value and the two-sigma range of an angle forthe normal population may be depicted in the image with three additional 
lines. Referencing may be either centrally, or with respect to one of the constituent measurement objects (that is with 
respect to one of the points in a point pair for distance measurements, or with respect to one of the lines in a line pair 
for angle measurements). This representation enables the clinician to graphically assess the actual measurement 
value in relationship to the average population. To avoid confusing superposition of too many lines in the radiographic 
image, all graphical mark-up pertaining to a particular measurand may be grouped in a separate layer, the display of 
which may be switched on and off on user's request. 

• Reporting and saving the results of a measurement session 

[01 24} Analogous to graphical and numerical input to a measurement session, both numerical and graphical results 
may be reported and saved. Measurement results and their applicable normal values are stored in the patient dossier 
in a suitable data format (such as XML). Furthermore, the graphic representation of a measurement session is saved 
as well in a suitable graphic format (such as SVG), and it is linked with the image, so that it may be redrawn onto the 
image at future referral. • 

Alternatively, the actualised graphics may be stored as a template, to be used in new measurement sessions during 
follow-up of a patient for example. 

General system overview 

{0125} A genera) overview of the measurement system is shown in fig. 1 . From a software functionality viewpoint, 
the system comprises atf engine and a user interface. • • 

(01 26} Theengine (E) guides the measurements that are performed on the retrieved (3) and displayed (4) radiological 
image by performing a generic execution of a measurement template. 

[0127] It loads andactivates a selected measurement template from a stored database (1) by (a) loading the meas- 
urement dependency graph Ihtomembry, (b) displaying (18) the graphical-part of the measurement template, (c) reg- 
istering all available methods specified by each measurement entity or object (both fextemar and Internal' methods 
associated with the graphical parr and the^gr^h hods part fespectryel)r>. 

It starts execution of the measurements at the topmost measurement entity as specified by.the measurement template 
pr as selected by the user, and traverses the measurement dependency graph whilechainlng the nodes along its path. 
When aleayenodeasioaaied with auserKiefmed measurement 

(5) of the point in the image. To this purpose, the online (a) passes contralto the user interface when the point mapping ' 
is performed bythe user who drags (6) the measurement objects onto an imageposifibh in accordance with the template 
after which the point's position is passed (7) to the measurement calculation module, or (b) initiates appropriate region 
of interest processing (e.g. contrast enhancement, geometric manipulation suc^a^-zobm , 
to enhance the desired anatomical landmarks), or to) performs the deformable contour-based segmentation for active 
measurement templates, in which case the position of the measurement points is established in spatial relationship to 
said contours and which positions are passed (7) to the measurement calculation module The initial position of the"' 
deformable contours may be determined using assail numberof mapped anchorpoirits (6). 
When art Internal node of the.graph fires (which is either a measurement object node or a measurement operator 
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node), the measurement computation module (A) computes the values of parameters of the notfe"baaed-an the values 
of the depending nodes. For example, the computation of the intersection point of two lines falls under this category. 
When the node is a measurement operator, the value of the measurement is computed and the measurement result 
is depicted graphically onto the radiological image (8). 

In a specific embodiment the engine further controls selection (9) and addressing (1 0) of normative values in a normative 
values table by a selection module (B), which normative values table is retrieved from a normal values repository (2). 
Both normative value (11) and actual measurement (12) are supplied to an abnormality judgement device (C), which 
calculates the presence ancTdegree of abnormality of the measurement value. The deviation from normality may be 
depicted graphically in the radiological image (13). The measurement results (14) and theirdegree of abnormality (15) 
is updated in a spreadsheet and clinical database. All graphical results of the measurement session are saved along 
with the image (16). 

[0128] The user interface (D) basically consists of a database, a controller and a graphics part. 
All graphical objects to be drawn in the image such as measurement objects (e.g. points, lines, circles...) and graphical 
measurement mark-up (e.g. distance and angle double-side arrow line segments) are generated and stored in a da- 
tabase created and maintained at run-time. 

The controller reacts to events such as the user pressing a mouse button or running a dialog box, starting a command, 
addition of an object to the database, notification and modification of objects that other objects onto which they are 
based have changed etc. 

The user interface's responsibility Is to render all measurement entities and measurement results. 

In particular it renders the measurement template in the template window (18), and changes appearance of entities in 

it in the course of the measurement procedure e.g. to hint the user to the next point to map. 

It generates or updates the graphical representation of all objects currently in the database, and superimposes it onto 
the display of the image in the image window (4). 

[01291 Itforther renders and updates the measurement values in the measurement values window (17). 

In a specific embodiment, normative values are displayed (20) in the measurement values window, along with the result 

of comparison of measurement values with their respective normative values (19). 

[01301 From the software architecture viewpoint, the system is organised into layers. * 
the system software layer, the bottom layer, contains components such as operating systems, databases, interfaces 
to-specific hardware and so on. , ; 

The middle-ware layer contains components such as GUI builders, interfaces to database management systems, plat- 
form-independent operating system services, and ActiveX-components such as spreadsheets and diagram ediiors. 
The business-specific layer, contains business specific components used in several applications. For example, a library . 
of calculations on geometric objects falls into this layer. 
The application layer, the top layer, contains the application specific services. 



Claims 

1 . A method of performing geometric measurements on a radiological image comprising the steps of 
■ - acquiring a digital Image representation of said radiological Image, 

- . applying said digital image representation to a display device connected to a computer and displaying said 
image, 

- activating a measurement scheme that has been stored in said computer, said measurement scheme con> 
prising a graphical part representing a relation between measurement entitles and anatomy present in said 
radiological image and an internal part representing functional dependencies between said measurement en- 
tities; said graphical part and said internal part being bi-directionatly linked, 

performing measurements on said displayed image, said measurements being guided by the activated meas- 
urement scheme. 

2. A method according to claim 1 wherein said measurements are performed by displaying said graphical part of said 
measurement scheme and by mapping measurement objects of the displayed graphical part onto a corresponding 
location in said displayed radiological image and by entering the corresponding location Into the activated meas-. 
urement scheme. • • ' • • ■ '.'- ■ 



A method according to claim' 1 wherein said measurement scheme is selected from a measurement scheme re- " 
pository comprising a number of measurement schemes adapted to different examination types and/or different 
radiological exposure conditions. 
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4. A method according to claim 1 wherein said measurement scheme is a passive measurefn^n^aeheme. 

5. A method according to claim 2 wherein constituent points of measurement objects are placed manually in said 
image. 

G. A method according to claim 2 wherein constituent points of measurement objects are placed by the steps of 
magnifying a region of interest surrounding the location of said constituent point in said image and/or in said meas- 
urement scheme and by manually placing said points in said image. 

7. A method according to claim 2 wherein values of measurands are computed by applying measurement operators 
onto geometric locations of measurement objects. 

8. A method according to claim 1 wherein said measurement scheme is an active measurement scheme. 

9. A method according to claim 1 wherein said graphical part is displayed as an overlay on top of the body part to be 
examined on said displayed image. 

10. A method according to daim 1 wherein said measurement scheme comprises normative values for measured 
entities and wherein upon activation of said measurement scheme said normative values are retrieved and com- 
pared with corresponding measured entities so as to assess a degree of abnormality. 

11. A method according to claim 1 wherein calibration factors are calculated for measurements based upon the image 
of a calibration device which is exposed simultaneously and under identical exposure conditions as the object to 
which said radiological image pertains, and wherein the measurements performed on said displayed image are 
corrected by means of said calibration factors. 

12. A method according to claim t wherein uncertainty values are calculated by error propagation of measurements 
and calculations guided by said measurement scheme 

-1-3. A method according to any of the claims 10, 11 or 12 wherein at least one of said calibrated values, uncertainty 
values, normative values and degrees of abnormality are visualised. 

1 4. A computer program product adapted to carry out the steps of any of claims 1 to 1 3 when run on a computer. 

1 5. A computer readable carrier medium comprising computer executable program code adapted to carry out the steps 
of any of claims 1 to 13. 
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